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Abstract. We consider the Cauchy problem for 2-D incompressible isotropic elastodynamics. 
Standard energy methods yield local solutions on a time interval [0, T/e], for initial data of the 
form cUq, where T depends only on some Sobolev norm of Uq. We show that for such data there 
exists a unique solution on a time interval [0, expT/e], provided that e is sufficiently small. This 
is achieved by careful consideration of the structure of the nonlinearity. The incompressible 
elasticity equation is inherently linearly degenerate in the isotropic case; in other words, the 
equation satisfies a null condition. This is essential for time decay estimates. The pressure, which 
arises as a Lagrange multiplier to enforce the incompressibility constraint, is estimated in a novel 
way as a nonlocal nonlinear term with null structure. The proof employs the generalized energy 
method of Klainerman, enhanced by weighted L 2 estimates and the ghost weight introduced by 
Alinhac. 



1. Introduction 

The long-time behavior of elastic waves for isotropic incompressible materials is studied in 2- 
D. The equations of incompressible elastodynamics display a linear degeneracy in the isotropic 
case; i.e., the equation inherently satisfies a null condition. By taking the advantage of this 
structure, we prove that the 2-D incompressible isotropic nonlinear elastic system is almost 
globally well-posed for small initial data. More precisely, we prove that for initial data of 
the form eU , there exists a unique solution for a time interval [0, exp(T(U )/e)], where T(U ) 
depends only on a certain weighted Sobolev norm of the Uq. 

To place our result in context, we review a few highlights from the existence theory of non- 
linear wave equations and elastodynamics. The initial value problem for small solutions of 3-D 
quasilinear wave equations with quadratic nonlinearities is almost globally well-posed [11], and 
in general this is sharp [8j. If, in addition, the nonlinearity satisfies the null condition, global 
existence was shown independently in [4J using conformal compactification and in [13] using the 
generalized energy method. The generalized energy method of Klainerman can be adapted to 
prove almost global existence for certain nonrelativistic systems of 3-D nonlinear wave equations, 
using scaling invariance to get weighted L 2 -estimates, as was first done in [13]. This approach 
was subsequently developed to obtain global existence under the null condition in [32] , see also 
[33] for a different method. A unified treatment for obtaining weighted L 2 estimates for certain 
hyperbolic systems appeared in [30] . 
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The existence question is more delicate in 2-D because, even with the null condition, quadratic 
nonlinearities have critical decay. A series of articles considered the case of cubically nonlinear 
equations satisfying the null condition, see for example [6j [T21 [21]. Alinhac [2] was the first 
to establish global existence for null bilinear forms. His argument combines vector fields with 
what he called the ghost energy method, but crucially it also relies on finite propagation speed 
through a certain Hardy-type inequality. 

The long-time existence theory for isotropic elastodynamics largely follows the paradigm of 
nonlinear wave equations. Almost global existence of small displacement solutions for 3-D 
compressible elastodynamics was first shown in [10], and counter-examples to global existence 
appeared in [9j [33]. The almost global existence proof was simplified considerably in [15] by 
enhancing the vector field approach with a new weighted L 2 estimate based on scaling invari- 
ance to compensate for the absence of Lorentz invariance. In [27], it was first noticed that 
there was a null structure compatible with the system of isotropic elastodynamics which can be 
used to establish global existence of small displacement solutions in 3-D. More comprehensive 
versions of this result appeared in [28, TJ. Compressible isotropic elastodynamics (in free space) 
is characterized by two wave families: fast pressure waves and slower shear waves. The afore- 
mentioned null condition limits the self-interaction of pressure waves. The equations possess an 
inherent null structure for shear waves. Thus, in the incompressible case where pressure waves 
are not present, it is plausible to expect global existence of small displacement solutions in 3-D 
without an additional null condition assumption, although the absence of finite propagation 
speed presents an obstacle. Nevertheless, this intuition was confirmed in [23 EI] . 

Our method for proving almost global existence for incompressible isotropic elastodynamics 
in 2-D is based on the ideas of [3T], with two new ingredients: the treatment of the pressure 
term and the use of the ghost weight. The equations are written as a first-order system with 
constraints whose unknowns are the material space-time gradient of the deformation, expressed 
in spatial coordinates (Section [2]). The advantage of this choice is that the Lagrange multiplier 
which enforces the incompressibility constraint appears as a pressure gradient. 

The generalized energy method forms the backbone of the argument. This enables control in 
L 2 (M?) of derivatives of the solution with respect to the vector fields which arise as infinitesimal 
generators corresponding to the fundamental invariance of the equations under translation, 
rotation, and scaling. This invariance gives rise to the basic commutation identities for the vector 
fields (Section HJ. The rotational vector field yields |x| -1 / 2 spatial decay, by means of Sobolev- 
type inequalities (Section [3]), which is weaker than the fully Lorentz invariant case. Using scaling 
invariance, the spatial decay is improved to t -1 / 2 time decay in L°°(IR 2 ) by means of a series of 
weighted estimates for the gradient of the solution which follow from an algebraic manipulation 
of the equations and the constraints. This algebraic procedure (which is an implementation of 
the method of [3D]) allows control of the solution gradient in L 2 (M. 2 ) with the weight (t — 
and it also gives t -1 decay in L 2 (IR 2 ) of certain special quantities (Sections El EJ. For small 
energy solutions, the weighted estimates are closed by a bootstrapping argument (Section [7J). 
The pressure gradient is estimated as a solution of a nonlinear Poissson equation (Section [5]). 
The nonlinearities of this elliptic equation have a null structure. Thus, the pressure term is 
essentially treated in a novel way as a nonlocal null form. This absence of finite propagation 
speed prevents us from achieving global existence. Energy estimates are performed using the 
ghost weight of [2] (Section E]) which provides a convenient solution to the technical problem 
that the weighted estimates hold only for the solution gradient. 

Finally, the general isotropic case is easily treated by our method since it can be regarded as 
a higher order nonlinear correction to the Hookean case (Section [TO]) . 
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Before ending the introduction, let us mention some related works on viscoelasticity, where 
there is viscosity in the momentum equation. In 2-D, the global well-posedness with near 
equilibrium states is due to [22] (see also [20]), and in 3-D, [3] and [19], independently. The 
initial boundary value problem is considered in [23] , the compressible case can be found in [261 E] • 
For more results near equilibrium, readers are referred to [3 I2U ESI E5J EE] ■ In [IE] a class of 
large solutions in two space dimensions are established via the strain-rotation decomposition 
(which is based on earlier results in [17] and [IB])- In all of these works, the initial data is 
restricted by the viscosity parameter. The work [13J was the first to establish global existence 
for 3-D incompressible viscoelastic materials uniformly in the viscosity parameter. 

2. Preliminaries and Main Results 

Classically, the motion of an elastic body is described as a second-order evolution equation in 
Lagrangian coordinates. In the incompressible case, the equations are more conveniently written 
as a first-order system with constraints in Eulerian coordinates. We start with a time-dependent 
family of orientation-preserving diffeomorphisms x(t, ■), < t < T. Material points X in the 
reference configuration are deformed to the spatial positions x(t, X) at time t. Derivatives 
with respect to the spatial coordinates will be written as (<9t,V) = d. Let X(t,x) be the 
corresponding reference map: X(t, •) is the inverse of x(t, ■). 

Lemma 2.1. Given a family of deformations x(t, X) , define the velocity, deformation gradient, 
and displacement gradient as follows: 

v{t,x) = ^^l , F(t,x) = ^p- , G(t,x) = F(t,x)-I. 

dt x=x(t,x) OX x=x(t,x) 

Then for < t < T, we have 

(2.1) djGik — dkGij = G m kd m Gij — G m jd m Gik, i,j, k e {1, 2, • • • , n}, 
and ifx(t,X) is incompressible, that is det F(t,x) = 1, then 

(2.2) V • v = dm = and (V • F T )j = d l F ij = 0. 

Proof. The first statement expresses the commutativity of material derivates DjDk = D^Dj in 
spatial coordinates. The second statement follows from Nanson's formula. Details are given in 

\m. □ 



Here and in what follows, we use the summation convention over repeated indices. The 
identities (12.11) and (12.21) will be used repeatedly in the sequel. 

To best illustrate our methods and ideas, we shall first consider the equations of motion for 
incompressible Hookean elasticity, which corresponds to the Hookean strain energy function 
W(F) = \\F\ 2 . 

{d t v + v ■ Vv + Vp = V ■ (FF T ), 
d t F + v-VF = VvF, 
V • v = 0, V ■ F T = 0. 

As will be seen in Section [10] where the case of general energy function is discussed, there is no 
essential loss of generality in considering this simplest case. Since the 3-D case has been treated 
in [291 EI] 5 below we will focus on 2-D. 
Denote the rotation operator by 

f2 = x 2 di — xi&2 = x x ■ V, 
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and the scaling operator by 

S = td t + Xidi + x 2 d 2 = td t + rd r , S = rd r . 
We shall frequently use the decomposition 

(2.4) V = (x/r)d r + (x ± /r 2 )n. 
Let r be any of the following differential operators 

(2.5) TE{d t ,d 1 ,d 2 ,Q,S}. 
As in previous works, we define 

£lf = V scalar /, 

Qv = Qv + (e 2 g) ei - e 1 <g) e 2 )v, V v G K 2 , 

hG = nG+ [(e 2 ® ei - ei ® e 2 ), G], V G G K 2 ® M 2 , 

where [A, 5] = AI? — denotes the standard Lie bracket product. 
Define the generalized energy by 

(2.6) E k (t) = J2 (lir a ^,-)lli 2 + l|r a G(t,-)||i 2 )- 

\a\<k 

We also define the weighted energy norm 

(2.7) X fc (t) = (ll(t-r)Vr%|| L2 + ||(t-r)Vr a G|| L2 ), 

|a|<fc-l 

in which we denote (a) = yl + a 2 . 

To describe the space of the initial data, we introduce 

A = {V,S ,Q}, 

and 

H k A = {f:J2 ||A a /IU» < oo}. 

\a\<k 

Then we define 

iZ*(T) = : [0,T) M 2 x (]R 2 <8> M 2 ) : (u, G) G n* =0 G^[0, T); 

with the norm 

sup E k (t) 1 ' 2 . 

0<t<T 

The main result of this paper is 

Theorem 2.2. Let (fo,Go) G u>rf/i > 5. Suppose that (vq,Gq) satisfy the constraints 
$11$, <0) and\\(v ,G )\\ K <e. 

There exist two positive constants Cq and €q which depend only on k such that, if e < eq, 
then the system of incompressible Hookean elastodynamics (12.31) with initial data (vo,Fq) = 
(vq, I + G ) has a unique solution (v,F) = (v,I + G), with (v,G) G H^(T), T > exp(G /e) and 
Ek(t) < 2e 2 , for <t <T. 
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3. L°° Decay Estimate 

In this section we derive several weighted L°° — L 2 Sobolev imbedding inequalities. These 
will be useful in proving decay of solutions in L°°. 

Lemma 3.1. For all radial functions f G H 1 (R 2 ), A = 1,2, there holds 
(3-1) Af(r)\ 2 <\\r^d r f\\ 2 L2 + \\f\\ 2 L2 

(3-2) r(t - r)"\f(r)\ 2 <\\(t- r)d r f\\ 2 L , + \\(t - r)*" 1 /^, 

provided that the right-hand side is finite. 

Proof. It suffices to show the lemma for / G Cq(R 2 ) and then use a completion argument. 
The first inequality is shown as follows 



poo 

r X \f(r)\ 2 = -r X d p [\f\p)\ 2 ]dp 

J r 

POO 

= -r x / 2f(p)d r f(p)dp 

J r 



oc 



< / P A - 1 |/(P)II^/(P)|P^ 
J r 

The calculation for the other inequality is similar. 

/CO 
d p [{t-p)'\f{ P )\ 2 ]dp 

/CO 
[(t- p ) x \t{p)\\d r f{p)\ + (t-p)^\f{p)\ 2 \dp 

/oo 
[(t-p) 2 \d r f{p)\ 2 +(t-p) 2 ^\f{p)\ 2 ]pdp 

= W-r)d r f\\ 2 L2 + \\(t-r)^f\\ 2 L2 . 

The next lemma is just the Sobolev imbedding theorem H 1 >■ L°° on the circle S 1 . 
Lemma 3.2. For each x = ru G IR 2 and f(ru>) G i/ 1 (S 1 ), there holds 



□ 



\f(rco)\ 2 <J2 [ \Wf{ru)\ 2 do. 



Proof. Write u = (cos 9, sin 9). We have 

| / fa) | 2 = | / fa) | 2 C os 2 9 + | / (ru) \ 2 sin 2 9 



r d 

= / — \\f(r cos <f>, r sin0)| 2 cos 
J* d<p L 

r e d 

+ / — \\f(r cos0, r sin 0)| 2 sin 
Jo d(p L 

[Iq/MH/mi + i/mi 2 ]^, 



< 

from which the result follows by the Cauchy-Schwarz inequality. □ 
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These two results combine to yield 
Lemma 3.3. For all f e H 2 (R 2 ), A = 1,2, there holds 
(3-3) r"\f(x)\ 2 < J2 W^WfWh + W^fWh 

a=0,l 

(3-4) r(t - r)"\f(x)\ 2 < £ [\\(t - r)d r n a ff L , + \\(t- rJ^OVIIl 

a=0,l 

provided the right-hand side is finite. 

Proof. First apply Lemma [3T21 to the left-hand side, then apply Lemma I3TT1 to the integrand. □ 

Remark. The estimates of Lemma 13.31 hold in the vector- and matrix- valued cases using Q in 
place of Q. This can be seen by applying Lemma 13.31 component-wise and the using the fact 
that, for example, \Qv\ < \Qv\ + \v\, for vectors v. 

Remark. The case A = 2 in (13. 4p will be applied only to derivatives V/. 
Lemma 3.4. For all f e H 2 (R 2 ), there holds 

(t)\\f\\L~ (r < m) < Yl \\(t-r)ff*f\\z>, 

\a\<2 

provided the right-hand side is finite. 

Proof. Let <p e C%°, satisfy <p(s) = 1 for s < 1, <p(s) = for s > 3/2. Note that (t) < (t - r) on 
supp (f(r/(t/2)). Thus, for |x| < (t/2), we have by Sobolev imbedding H 2 (R 2 ) ■=— >■ L°°(1R 2 ) that 

(t)\f(x)\ = (tMr/(t/2))\f(x)\ 

<(t)y(r/(t/2))f\\ H2 

< (t) E ll^/H^(,<|(V2» 

|o|<2 

< W-r)d a f\\ L >. 

\a\<2 

□ 

Remark. This result will only be applied to derivatives V/. 

4. Commutation 

For any multi- index a, we have the following commutation properties when applying T a to 
the equation (12. 3p (see, for instance, [3TJ [13] for details). This will be essential in all of the 
subsequent estimations. 

'd t T a v - V ■ T a G = -Vr> + f a , 
(4.1) { d t T a G - VT a v = g a , 

V • T a v = 0, V • T a G T = 0, 



where 



(4.2) 



/«= E [-^■vn ? ; + v.(PGr^ T )], 

/3+7=a 

3a = E [-T p vVnG + VnvT p G], 

/3+7=a 
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From (|2.ip . we also have 

(4.3) V 1 ■ r Q c = /i Qj 
where 

(4.4) (/i a ) i= ]T [r^^r^-r^ATOa]. 

,3+7=0 

5. Bound for the Pressure Gradient 

The following Lemma shows that the pressure gradient may be treated as a nonlinear term. 
The first estimate appeared in [31 J . The second is a novel refinement which saves one derivative 
over the first bound and which allows us to exploit the null structure. This is essential in Section 
[9] when we estimate the ghost weighted energy. 

Lemma 5.1. Let (v, F) = (v, I+G), (v, G) E H^(T), solve the equation (12.31) and the constraint 
( 12. ID . Then we have 

(5.1) l|vr>|| L2 < ||/ Q || i2 

(5.2) l|vr>|| L2 < \\d j r p v i r'v j -d j rf > G ik rtG jk \\ v ,, 



p + 7 = a 

\p\ < h\ 



for all \a\ < k — 1. 



Proof. Taking the divergence of the first equation of (14.11) and then using the constraints given 
in third equation of ( 14. ip . we find 

AT a p = V • f a + V • (V • T a G) - d t V -T a v = V- f a . 

By (14.21) and the constraint equations in (14. ip . we have 

v-/ a = - Yl did j (r0v i r'v j -rf l G ik r'G jk ) 

P+"/=a 

d i d j (Y fi v i r'v j -T fi G ik r'G jk ) 

P + 7 = a 

\p\ < h\ 

did^v^v^V^G^G^) 

P + 7 = a 

\P\ > h\ 

= d i (d j Y fi v i r'v j -d j T fi G ik r'G jk ) 

P + 7 = a 

\p\ < h\ 



d j (T p v i d i rv j -r p G ik d i r r, G jk ). 



P + 7 = a 

\p\ > h\ 



The result now follows since 



vr> = a- x v(v • f a ) 

and since A _1 V <S> V is bounded in L 2 . 



□ 
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6. Estimates for Special Quantities, I 
Lemma 6.1. Suppose that (v,F) = (v,I + G), (v,G) G H k K (T), solves Q, flU]). Define 

L k = ^ [|r%| + 

|a|<fe 

and 

(6.1) N k = [t\M+t\9a\ + (t + r)\h a \+t\Vr a p\]. 

\a\<k-l 

Then for all \a\ < k — 1, 

(6.2) r\d r T a vu\ <L k 

(6.3) r\d r T a G T u\<L k 

(6.4) r|9 r r Q Gw - V • T a G\ < L k 

(6.5) r\d r T a Gio L \<L k + N k 

(6.6) (t±r)|Vr a t;± V-TO®w| <L fc + AT fc . 

Proof. By the decomposition (12 ,4p we have, 

(6.7) VT a v = d r T a v <g> to + -QT a v ® uA 

r 

Taking the trace of this identity yields 

V • T a v = d r T a v ■ lo + -QT a v ■ u L , 

r 

and so, by the divergence- free velocity constraint of (14.11) . we obtain (I6.2p . It also follows from 
flSTTj) that 

(6.8) r\VT a v - d r T a v ® w| < L fc , 

which will be used shortly in proving (16. 6p . 

Again by (12. 4ft . we may write for any matrix- valued function H 

(6.9) d r H = d r H I 

= d r H[u g) to + to 1 - (g) w 1 "] 
= d r Hco <g> oj + d r Hco 1 ' <g> w -1 

= [V • if - -OHu^] <g> w + [V x • # + iflifw] ® uA 
Multiplying both sides of (16 .9p times the vector w, we obtain 

<9 r Fu; — V-H— -QHu L . 

r 

Apply this to H = T a G T , and use the other divergence- free constraint from (14. ip . We thereby 
obtain 

(6.10) r\d r T a G T u\ = |Or Q GV| < L k , 
which is (16.31) . 

Next, apply (16. 9p to H = T a G and use the constraint (14. 3p : 

rV ■ T a G ® u - rd r T a G = QV^Gu 1 ® u - [rh a + ttT a Gcu] ® u ± . 
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From this it follows that 

r\V-T a G-d r T a Guj\ = r|[V • T a G ® u - d r T a G}u\ = \nr a Gu x \ < L k , 
proving (16. 4p . and also 

(6.11) r\V -T a G®uj-d r T a G\ < L k + N k . 

As an immediate consequence of this last inequality, we get 

r\a r r a Gu J -\ =r|[V ■r a G®u-d r T a G]u x \ <L k + N kl 

which proves (16.51) . We are now ready to prove (16.61) . 

Using the PDEs (14. ip and the definition S = tdt + rd r , we can write 

tVT a v + rd r T a G = ST a G - tg a 

tV • T a G + rd r T a v = ST a v - tf a + tVT a p. 

This is rearranged as follows: 

tVT a v + rV • T a G g> u = r[V ■ T a G ® w - <9 r r Q G] + ST Q G - £# a 

tV ■ T a G ® w + rVT a v = r[VT a v - d r T a v <g) w] + [Sr Q i; - t/ a + tvr°p] ® w. 

Notice that by (16.81) and (16. lip , the right-hand sides of these identities are bounded by L k + N k . 
Therefore, the bounds (16. 6p result from the combination of these two identities. □ 

Lemma 6.2. Let (v, G) G H\(T). Then for \oc\ < k — 2, we have 



il/2 

Proof. We shall make use of the fact that 



rT a v ■ lu\\ L oo + \\rr a G T u\\ L °o < E^ + . 2 . 



(6.12) n(v(x)-u) = (nv(x))'U>, n(G(x) T cu) = (QG(x)) T u, cu = x/r. 

By ([331) with A = 2 and (16321) . we have 
||rr Q va;|| L oo + ||rTO T u;|| L oc, 

< ^ [\\rd r Q a {T a v ■ u)\\ L 2 + \\rd r Q a {T a G T u)\\ L 2 



o=0,l 



+ \\n a (T a v ■ oo)\\ L 2 + \\n a (T a G T uj)\\ L 2} 
= [\\r{d r Q a T a v) ■ u\\ L 2 + \\rd r {Q a T a G T )u\\ L 2 

o=0,l 

+ ||(fi a r%)-u;|| L2 + ||(fi a r Q C7 T )a;|| L2 ]. 

The result now follows by (16.21) . (16.31) . □ 

7. Weighted L? Estimate 
In this section, we show that the weighted norm is controlled by the energy, for small solutions. 
Lemma 7.1. Suppose that (v, F) = (v, I + G), (v, G) G H k k (T), k > A, solves (|2T3]) . (l2TTj) . Then 

\\N k {t)\\ L 2<E k {t) + E k {t) l l 2 X k {t) l l 2 . 
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Proof. By (16. ip and Lemma 15.11 we have that 

\\N k {t)\\& < E Wir a p(OII^ + *ll/a(i)IU 2 + ^II^WIU 2 + ||(t + r)/i Q (t)|U 2 ] 

\a\ < k- 1 

< E [*ll/a(*)ll^+*IM*)llL a + ll(* + » , M*)IW 
|a| < k- 1 

In estimating these terms, we shall consider two regions: r < (t/2) and r > (t/2). 



Estimates of nonlinearities for r < (t/2). Examining definitions (14. ip . f)4.4p . we find that 



E WI./a|U 2 (r<<t/2» +t\\g a \\ L 2 (r<{t/2)) + \\{t + r)h 

\a\ < k- 1 

< E (*)ll(|r^i + |r^|)(|vr^| + |vTO|)|| i2(r < (t/2)) . 

/3 + 7 = a 
\a\ < k- 1 

To simplify the notation a bit, we shall write 

\(v k vj k G)\ = ^[|r Q ^| + |r a G|], 

|a|<fc 

(and similar). We make use of the fact that since /3 + 7 = a, \a\ < k — 1, and k > 4, either 
+ 2 < k or | 7 | + 3 < k. Thus, we have by (E3D 

(t)||(|r^| + |r^|)(|vr^| + |vTO|)|| i2(r < (t/2)) 

<<*> II (r k ~\ r k - 2 G) 11 (vr fc - V vr^G) \\^ {r < m) 

+ (t) || (T k v, T k G) \\ L i I) (VT k - 3 v, VT k - 3 G) \\ L °°(r< m) 

<\\ (r k - 2 v, r k - 2 G) \\ L ~ 11 (t - r) (vr fc -y vr^G) \\ LHr < m) 
+ \\(r k v,r k G)\\ L 2 \\(t-r)d\vr k -\vr k - 3 G)\\ L ~ (r < m) 

|A|<2 

<E k {ty'*x k {tyi*. 



Estimates of nonlinearities for r > (t/2). Using (12.41) . we replace all derivatives which 
occur in (14.21) . (14.41) by their radial and angular parts. We find that 



\f a \ + \9a\ + \h a \ < E (l r ^ ■ u \ + \^G T u\){\drT 1 v\ + \d r VG\) 

+ 7 = a 
\a\ < k- 1 

+ (\T p v\ + \r p G\)-(\nr r v\ + \nr f G\). 

/3 + 7 = a 
\a\ < fc - 1 
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Thus, we obtain 

I 7 - 1 ) [*ll/«||L2(r>(t/2» + t\\g<x\[L*{r>{t/2)) + \\(t + r)h 

a \ \L 2 (r>(t/2))\ 

|«|<fc— 1 

< ll r d r ^ • u \ + |r^ T o;|)(|a r r^| + \d r VG\)\\v {r > m) 



$ + 7 = a 
\a\ < fc-1 



+ E IKI r ^l + \^G\)(\nr^v\ + \n^G\)\\ LHr 



><*/2»- 



/3 + 7 = q 
|a| < jfe- 1 



We claim that all terms on the right-hand side of (17.11) can be bounded by E^. 

Consider the term in the first sum which has — a, 7 = 0. It can be estimated as follows: 

||r(|r a v ■ a;| + |r a G T a;|)(|^| + \d r G\)\\ L . {r > m) 
(7.2) < ( sup r 2 [ (\T Q v(ru) ■ co\ 2 + \T a G T (ru)u\ 2 )da ) 

\r>(t/2) Js 1 J 

\ 1/2 

x( / sup(\d r v(rcu)\ 2 + \d r G{ruj)\ 2 )rdr 

\J(t/2) §1 

Using ( 13. II) with A = 2 from Lemma 13.11 the first of these integrals is bounded by 
\\rd r T a v ■ u\\ L 2 + \\rd r T a G T Lo\\ L 2 + \\T a v\\ L 2 + ||r Q G|| L 2. 

Noticing the form of the first two terms, we can use 

flHZ|,fE3D to bound this by E^ 2 +1 < E]' 2 
The second integral in ( 17. 2ft is bounded as follows using Lemma [3.21 

£ m a drv\\ L 2 + \\tt a d r G\\ L .) < e\ /2 

a=0,l 



k 



This proves that the term in (I7.2p is bounded by E^, as claimed. 

For the remaining terms in the first sum in (17.11) . we have /3 7^ a. We estimate these as 
follows: 

(7.3) ||r(|r^ ■ u\ + \T^G T u\)(\d r T^v\ + \d r r>G\)\\v (r > m) 

< (\\rV p v ■ ouWlo. + ||rr /3 C7 T a;|| L oc)(||9 r r^|| L2 + \\d r T^G\\ L2 ). 

The terms in L°° are estimated by E,L +2 , using Lemma 16721 and thus, since < |a| — 1 < k — 2, 
the expression ( 17. 3ft is bounded by 

T-,1/2 771I/2 7-1 

Altogether, the first sum of terms in ( 17. lft is bounded by Ef.. 

The second group of terms on the right of ( 17.11) is also bounded by Ek, using the same 
interpolation strategy as in the case r < (t/2). 

□ 

Lemma 7.2. Suppose that (v, G) G H k K (T), solves (Q, (J2HJ- 

|| (t - r)Vr a G|| L2 < E k {tfl 2 + \\{t- r)V ■ r Q G|| L2 + \\{t + r)h a \\ L i. 
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Proof. For any IR 2 ® M 2 - valued function H, we have that 

\VH\ 2 - [|V • H\ 2 + |V X • H\ 2 } = -2[a 1 ( J Ff a 9 2 F i2 ) - ^(HadxH^)]. 
Thus, using integration by parts and Young's inequality, we get 

||(* - r)VH\\l 2 - [\\(t - r)V • H\\\ 2 + ||(f - r)V^ • H\\ 2 L2 ] 

= - [ 2[d 1 {H iX d 2 H i2 ) - d 2 (H il d 1 H i2 )]dx 
<l\\(t-r)VH\\ 2 L2 +C\\H\\l 2 , 

and so we obtain 

IK* - r)VH\\ 2 L2 <\\(t- r)V ■ H\& + \\(t- r)V x ■ ff|& + ||tf || 2 2 . 
The Lemma follows by applying this inequality to H = T a G and then using (14.31) . 



□ 



Theorem 7.3. Suppose that (v,F) = (v,I + G), (v, G) e H\{T), k > 4, solves (E3D, (EITJ. If 
E k (t) < 1, thenX k (t) < E k {t) 1 / 2 . 

Proof. Starting with definition (12 .7p and using the fact that (t — r) < 1 + \t — r\, we obtain from 
Lemma 17.21 

x%= E [IK*-^)vr^||i 2 + ||(t-r)vr a G||i 2 ] 

i«i<fc-i 

<E k + J2 [ll(*-r)vr^||i 2 + ||(t-r)vr-G||i 2 ] 



H<fc-i 



< 



E k + E [||(*-r)Vr^|ll 2 + ll(*-r)V-TO||| 2 + ||iV fc ||| 2 ]. 



\a\<k-l 



Since 



and 



vr°w = ^[vr a v + v • r a G ® u] + ^[vr°w - v ■ r Q c <g> u] 
v • r a G = -[vr a v + v ■ r a c ® w]w - i[vr Q v - v • r Q c <g> 



we see that 
\t-r\ [\VT a v\ + |V-r Q G|] 



<\t + r\ | Vr Q f + V • T a G <g> cj\ + \t - r\ \VT a v - V • T a G ® w|. 

It follows from f 1 6 . 6 f) that 

||(f - r)Vr^||| 2 + ||(* - r)V • r-C||| 2 < E k + ||iV fc || 2 . 

and thus, we obtain 

x fc <^ /2 + ||Ay L2 . 



Then applying Lemma [7. 11 we get 



x fc <^ /2 +E fc +£;y 2 x t i/2 
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Under the assumption that E^. < 1, we obtain 

X k < E l J 2 + E k < E l J 2 . 

□ 



8. Estimates for Special Quantities, II 

With the results of the previous section in hand, we can now complete the estimation of 
T a v + T a Gu and T Q Gw ± . 

Lemma 8.1. Let k > A, E k -C 1, u — x/\x\. Then we have 

(8.1) \\r(d r T a v + d r T a Guj)\\ L 2 + \\rd r Y a Guj^\\ L 2 < E l k /2 , \a\<k- 1, 

(8.2) ||r(r^ + r a L7w)|| L oo + ||rr Q C7a; ± || L o < E l k ' 2 , \a\ < k - 2. 
Proof. First, note that 

d r T a v + d r T a Guj = (VT a v + V ■ r Q C7 ® u)u + (d r T a Guj - V • r°C7). 

Therefore, by (16. 4p and (16.61) . we have 

r\d r T a v + d r T a Gio\ < (r + t)\VT a v + V • T a G g> u\ + r\d r Y a Gu - V • r Q G| < L k + N k . 

Combining this with (|6.5p . we get the estimate 

\\r(d r T a v + d r T a Gco)\\ L 2 + \\rd r Y a Gu L \\ L 2 < E^ 2 + \\N k \\ L 2. 

Estimate (18. ip now follows by Lemma 17.11 and Theorem 17.31 

To obtain the other estimate, we observe that (similar to (I6.12p ) 

(8.3) Q{T a Gu) = {QT a G)u, tt^Gu^ = (nr a G)u ± . 
By ([33D with A = 2, (jO]) . flHTTD . we have 

\\r(T a v + r a C7w)|| L oo + \\rT a Guj L \\ L o. 

< [\\r{d r tt a r a v + d r Q a T a Gu)\\ L 2 + \\Q a T a v + Q a T a Gu\\ L 2 

a=0,l 

+ \\rd r Q a r a Gu L \\ L 2 + \\£V t r a Gu ± \\v) < E l k ' 2 . 

□ 

9. Energy Estimate with a Ghost Weight 

Choose q — q(t — r), with q(o~) = f°(z)~ 2 dz so that q'(a) = (er)~ 2 and \q(o~)\ < |. Let 
|a| < k. Taking the inner product of the first and second equation in ( ]4.ip with e~ q T a v and 
e~ q T a G respectively, and then adding them up, we find 

J [e- q d t (\T a v\ 2 + |r Q Gf ) - 2e- q (T a v i d j T a G ij + T a G i &T a v i ))dx 

= -2 / e-r> • to + 2 / .-«</. • i-» + (9 j,r« G ,)^. 
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Integration by parts gives that 

(9.1) 4 I e- q (\T a v\ 2 + \T a G\ 2 )dx 
at J 

= "/ (T^ (|r °" + r ° Gw|2 + l r " &± l 2 )^ 

- J / e ^..vr^ + ,/.^.^ + ( fc )^ 

We emphasize that here we do not use integration by parts in the term involving pressure. We 
also point out that we cannot use the approach of Lemma 17.11 to estimate the nonlinear terms 
because we now have that \a\ < k rather than \a\ < k — 1 as we had earlier. 

Let us first treat the last term in (19. ip . Recall that f a and g a are given by (14. 2p . Since v and 
G T are divergence-free, we get 

(9.2) J e-" [f a ■ T a v + {g a ) i3 Y a G i3 ] dx 



P + 7 = a 
1 + <* 



P + 7 = a 



+ \J e-^[2r^r a G lfc G Jfe -^(|r^| 2 + |r Q G| 2 )]rfx. 

To estimate the last term in (19. 2p . we first compute that 
J e'ldj [2T a v % T a G %k G jk - v,{\T a v\ 2 + \T a G\ 2 )]dx 

= - J{t- r)- 2 e- q [2T a v i r a G ik G jt fu j - (\T a v\ 2 + |r Q G| 2 ) W ] dx. 

This is estimated by 

(|| (t - r)- 2 v -cu\\ L ~ + \\(t - r)- 2 G T u\\ L oo)E k (t). 



Now by Lemma I6.2[ we have 

|| (t - r)~ 2 v ■ u\\ L oc < \\{t- r)~ 2 v ■ w|| L oo( r <( t /2» + IK* - r )~ 2 ^ ' ^\\h°°{r>(t/2)) 

< (t)~ 2 |M|L°°(r<(</2» + \\V ■ v\\L°°(r>{t/2)) 

^ (*)~ 2 |MU°°(r<(</2» + (*}~ || rl ' ' w IU°°(r>(t/2» 

< (t)^E 1 ' 2 . 

A similar estimate holds for the term with G T u. We have shown that the last term in (19.21) is 
bounded by (t)' 1 ^ 2 . 
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Next, we are going to estimate the first and second terms on the right hand side of (19.21) . 
Since k > 5, it is always the case that |7| < k — 1 and either < k — 2 or I7I < k — 3. 
For r < (t/2), we use ( 13.4p and Theorem 17.31 to get 

(9.3) J2 I / e-*r a v i [d j r'Gi k r fi G jk -r f) v j d j r'v i ]dx 

„ 1 Jr<(t/2) 



(3 + 7 = a 
+ 



'r«t/2> 

/ e-«r a G ik [d^v^G.k - Y^v j d j r'G ik ] dx 

Jr<(t/2) 

< (t)- 1 ^^) 1 /^!! (t _ r) (vr fc -^, vr^ 1 ^) |U 2(r .« V2)) || (r fc ^, r fc - 2 c) ||^ 
+ ||(t)(vr fc - 3 w ,vr fc - 3 G)|| L oc (r<<t/2)) ||r fe (t;,c7)|| L2 



< {t)- l E k (tf\ 
For r > (t/2), we write using (12. 41) 



(9.4) 



V I / e- q Y a v i [d j r'G (k Y p G ik -Y p v i d i r'v i ]dx 

1 Jr>{t/2) 



/3 + "/ = a 
7/« 



+ 



< 



/ e-"Y a G ik [d^ Vi Y^G jk - Y^ Vj d^G lk ] dx 

Jr>(t/2) 

V I / e-^VilujjdrY^GikY^Gjk-Y^v^jdX'v^dx 

1 Jr>(t/2) 



/3 + 7 = a 
1 + a 



+ 
+ 



/ e~ q Y a G ik [ujdrVviT^Gjk - r^-u^r 7 ^] dx 

Jr>(t/2) 



fl f \ doc 



r>{t/2) 



V I / e- q ((Y a v,Y a G),B[(Y^v,Y^G),(d r Y^v,d r Y^G)})dx 

1 Jr>{t/2) 



(3 + 7 = a 
7 ^ a 



+ 



-fijj-j/ doc ■ 



r>{t/2) 



in which 



and 



B[(u x , Gi), (v 2j G 2 )] = (G 2 Gju - ( Vl ■ u)v 2 , v 2 ® G^u; - ( Vl ■ u)G 2 ), 



/3 + 7 = a 
7 ^ a 
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To analyze the structure of these terms when r > (t/2), we decompose vector/matrix pairs 
(v, G) as 

1 

(v,G)= ^iL>,G) 

k=-l 

Utiv, G) = \{{v + Gcj), (v + Gcj) <g> w)) 
n_i(u, G) = - Gu;), -(« - Gw) (8) w)) 
n (t>,G) = (0,G?w- L ®w L ). 
We have the following cancellations 

B[n fc (v 1J G 1 ) J n fc KG 2 )] = o, fe = ±i 

B[EIo(vi,Gi),n fc (v 2j G a )]=0, fc = ±l 

5[n fe K, Gi), n (t; 2 , G 2 )] =0, k = ±1. 

In particular, there is no self- interaction of the "bad" quantity II_i. Thus, we have 

B[( Vl , Gi), (v 2 , G 2 )] = B^iv!, Gi), n_!(t; 2> G 2 )] 

+ Gi),IIi(v2, G 2 )] + B[U ( Vl , Gi), n (^ 2 , G 2 )]. 

The conclusion of this discussion is that the terms in (19 Ah are bounded by 



+ 7 = a 



|(r Q v,r a G)| |r^ + r^Gw| |(<9 r r 7 v,«9 r r 7 G)| 



r>(t/2> 



+ |(r^,r^G)| |9 r r 7 w + 9 r r 7 Gu;| 
+ -|((r^,r^G)| \(nr^v,hr^G)\ 

r 

<C\\(T k v,T k G)\\ L . 



dx 



x 



(^-Ir^-^ + r^G^iUoo^^j^T^,^-^)!!^ 



-1/2 I 



^r 1 ' 2 ^ - r^drT^v.drT^G)^^) 



(t - r) 

+ (trKr^-y r fe - 2 G) || ioo ||r|a r r fc - 1 i; + air^Gwi + r|d r r fe ~ 1 Gu; ± ||U 2{r>(t/2>) 
+ (t)-l(ryr fc G)|| L2 ||r|«9^ 



By (18. 2p . we have that 
By (13. 4p . we obtain 

1/a r^ + r fc Ga; 



(t-r) 



L 2 



|r 1 /2(t_ r )( 9r r fe - 3 ? ;,a r r fc - 3 G)|| I/ o 0(r>(V2>) 



<<*> 



-1/2 



r fc w + r fe Gw 



E 



1/2 



L 2 
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By the conventional Sobolev imbedding, (|8.ip . and f 18 . 2 j) we have 



r fe - 2 V) r fe - 2 G)|| L o ||r|a r r fc - 1 ? ; + a r r /c - 1 Ga;| + r|a r r ,c - 1 G'a; x ||| L 2 < £ 



fc-l. 



and 



||(r fc v,r fc G)|| La ||r|a r r fc - 3 v + a P r fc - 3 Ga;| + r|a r r fc - 3 Ga; ± ||| L o <E k . 

It follows, therefore, that the integral (19 .2p is bounded by 



(t-r) 



L 2 



for an arbitrarily small > 0. 

It remains to treat the pressure term in (19.11) . However, thanks to (I5.2p . this term is handled 
exactly as the preceding ones. 

Finally, we gather our estimates for (19. ip to get 



m + 



T k v + T k Gu 



(t-r) 



L 2 



T k v + T k Gu 



(t-r) 



L 2 



+ C,(t)^E 3 k /2 , 



with 



E k {t) = I e-\\T a v\ 2 + \T a G\ 2 )dx. 

\a\<k 



Noting that E k ~ E k , we obtain, for fi small, 

E' k (t) < C7 (t)- 1 ^(t) 3/2 , < t < T. 
This implies that E k {t) remains bounded by 2e 2 on a time interval of order T ~ exp(Go/e). 

10. General Isotropic Elastodynamics 
For general isotropic elastodynamics, the energy functional has the form W = W(F) with 
(10.1) W(F) = W(QF) = W(FQ) 

for all rotation matrices: Q = Q T , det Q — 1. The first relation is due to frame indifference, while 
the second one expresses the isotropy of materials. This implies that W depends on F through 
the principal invariants of FF T , namely tr FF T and det FF T in 2-D. Setting r = | tr FF T 
and 8 = det F = (det FF T Y^ 2 , we may assume that W(F) = W(t, S), for some smooth function 



W : R + x 
Since 



— > 



- = F and W = SF-', 



the Piola-Kirchhoff stress has the form 

dW(F) 



S(F) 



dF 



W t {t,5) F + W 5 (t,S)5F- t . 



We assume that the reference configuration is stress free, S(I) = 0, so that 

(10.2) W- T (1,1)+W*(1,1) = 0. 

The Cauchy stress tensor is 

T(F) = 5- 1 S{F)F T = 5~ l W T (r, 5)FF T + W s {t, 5)1, 

and the term V • FF T in (12.31) is replaced by V • T(F). Let us now proceed to examine this 
term. 
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Write 

T(F) =W T (1, 1)FF T 

+ {6- 1 W T (r,5) -W T (IA)}[FF T - 1} 
+ [8- l W T {r, 5) - W T (1, 1) + W s (r, S)]I 

0=1 

Assume that 

(10.3) W T (1,1)>0. 

Then Ti(F) gives rise to a Hookean term. Notice that assumption fl 1 . 3 1) rules out the hy- 
drodynamical case W T = 0. The principal invariants can be expanded about the identity as 
follows: 

(10.4) r = -tr FF T = ^tr (/ + G)(I + G T ) = 1 + tr G + -tr GG T 
and 

(10.5) 5 = det F = det(J + G) = 1 + tr G + det G. 

In the case of incompressible motion, we have 5 = 1, so from ( 110. 5p . we get 

tr G + det G = 0, 

and hence from fjl0.4j) 

t - 1 = - tr GG T - det G = 0(\G\ 2 ). 
Thus, we see that for |G| 1 

T 2 {F) = 0[{r-l)\G\] = 0{\G\\ 

which produces nonlinearities which are of cubic order or higher. Finally, T^{F) leads to a 
gradient term which can be included in the pressure. The conclusion is that the general incom- 
pressible isotropic case differs from the Hookean case by a nonlinear perturbation which is cubic 
in the displacement gradient G. Such terms present no further obstacles in the proof of almost 
global existence in 2-D, in particular, they have the requisite symmetry properties for energy 
calculations, see [3T] . 

Theorem 10.1. Let (t> ,G ) € H*, with k > 5. Suppose that (v ,G ) satisfy the constraints 
dSD, tfOD and ||(« ,Go)||H*<e. 

Assume that the smooth strain energy function W(F) is isotropic, frame indifferent, and 
satisfies pU2!) . ffTOTBl . 

There exist two positive constants Cq and €q which depend only on k such that, if e < eo, then 
the system of incompressible isotropic elastodynamics 

{d t v + v ■ Vv + Vp = V • T(F), 
d t F + v-VF = VvF, 
V • v = 0, V • F T = 0. 

with initial data (vq, Fq) = (vq, I + Go) has a unique solution (v, F) = (v,I + G), with (v , G) G 
Hl(T), T > exp(C /e), and E k (t) < 2e 2 , for0<t<T. 
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